Tunnelings as Catastrophes 



C. A. A. de CarvalhcQ 

Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
Cx. Postal 68528, CEP 21945-970, Rio de Janeiro, RJ, Brasil 

R. M. Cavalcantfl 

Departamento de Fisica, Pontificia Universidade Catolica do Rio de Janeiro, 
Cx. Postal 38071, CEP 22452-970, Rio de Janeiro, RJ, Brasil 



in 

On 
On 

> 

o 
m 

(N 
> 

in 
o 
o 



in 

ON 



X 



We use path-integrals to derive a general expression for 
the semiclassical approximation to the partition function of 
a one- dimensional quantum-mechanical system. Our expres- 
sion depends solely on ordinary integrals which involve the 
potential. For high temperatures, the semiclassical expres- 
sion is dominated by single closed paths. As we lower the 
temperature, new closed paths appear, including tunneling 
paths. The transition from single to multiple-path regime 
corresponds to well-defined catastrophes. Tunneling sets in 
whenever they occur. (Our formula fully accounts for this 
feature.) 
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We may use path-integrals in the description of Quan- 
tum Statistical Mechanics. The partition function for 
a one-dimensional quantum-mechanical system, for in- 
stance, can be expressed as 



Z(J3) 
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For an arbitrary potential, V(x), this integral has been 
approximated by means of perturbativc jlj|^] , variational 
jl],U and numerical ||] techniques. Here, we shall concen- 
trate on the semiclassical approximation to the integral 
in order to: i) derive a general formula for Z sc (f3) that 
does not require a detailed knowledge of the classical Q] 
motion; ii) discuss the onset of tunneling and relate it to 
the study of singularities in Catastrophe Theory. 
The semiclassical evaluation of (|la| ) yields: 

Z SC (P) = / dx Y, e- SiM)/R AJ 1/2 (xo,P), (2) 

where Sj(xo,f3) denotes the action and Aj(xo,(3) repre- 
sents the determinant of the fluctuation operator, 



(3) 



both calculated at the j-th classical trajectory, Xj(r), 
satisfying the boundary conditions, Xj(0) = Xj({3H) = xq; 



N(xo, 0) is the number of classical trajectories which are 
minima of the action functional. 

The action has a simple expression in terms of the turn- 
ing points of the classical trajectory: 



Sj(x ,/3) = (3hV(x j ± )±2 / dxMv(x,x 3 ±) 



L 



+ 2n / dx Mv(x,x 3 ± ) 



(4) 



Here, v{x,y) = y/(2/M)[V(x) - V(y)) and x\ (xt) are 
turning points to the right (left) of xq . The first term in 
([|) corresponds to the high-temperature limit of Z((3), 
where the classical paths collapse to a point, i.e., — > 
xq. For motion in regions where the inverted potential 
is unbounded (hereafter called unbounded motion), n — 
0. For periodic motion, n counts the number of periods 
and the second term is absent. For bounded aperiodic 
motion, all terms are present. The last two terms will be 
negligible for potentials which vary little over a thermal 
wavelength, A = Tiy/ (3/M. However, by decreasing the 
temperature, quantum effects become important. 

For trajectories having a single turning point (n = 0), 
x± are given implicitly by 



t(x ,x 1 ± ) 



dx 



Ixo v(x,X±) 

and the fluctuation determinant by 
4Trh[V(xi) - V(x )} 



Aj(x o ,0)=±- 



MF'(4) 



dT(x ,y) 
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The formula above is valid only for trajectories with a 
single turning point. However, this is not really a restric- 
tion, since trajectories with two or more turning points 
are naturally excluded from our calculations, as it will 
become clear in the sequel. 

Formulae ([|) and (|^) do not require full knowledge of 
the classical trajectory, as the dependence on xj (t) comes 
only through the turning point. To see how this can 
simplify the evaluation of Z sc (/3), let us take the harmonic 
oscillator, V(x) = ^Maj 2 x 2 , as an example. In this case, 
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given Xo and (3 there is only one trajectory, with a single 
turning point, given by x+{x-) — xq/ cosh(/3hu> / 2) for 
x < (> 0). Si and Ai can also be readily calculated, 
the final result being 



Z sc (l3) 



dx e -( M ^o/ s ) tanh (/3^/ 2 ) 



Mlo 



2nh smh(/3huj) 



(7) 



which in this case is exact. 

A more interesting situation occurs in the case of the 



anharmonic oscillator, V(x) — \(x 



Jl\2 



For x 2 > a 2 



only single paths with single turning points exist for fixed 
xq and /3. However, there is also a region, x 2 < a 2 , 
where the classical motion is bounded and a much richer 
structure exists, one in which more than one classical 
path may exist for given values of xq and (3. 

In a region of bounded classical motion (a well in —V), 
such as x 2 < a 2 for the anharmonic oscillator, the number 
of classical trajectories changes as the temperature drops. 
If < /3 < ir/h,uj m (where ui m = yj— V"(x m )/M and x m 
is a local minimum of — V), for every xq in this region 
there is only one closed path, with a single turning point, 
satisfying the classical equations of motion. For xq < 
x m (> x m ) this path goes to the left (right) and returns 
to xo. For xo — x m , it sits still at the bottom of the well. 
It is this single-path regime which goes smoothly into the 
high-temperature limit. 

For (3 — ir/huj m , the solution that sits still at x m be- 
comes unstable. Its fluctuation operator is that of a har- 



monic oscillator with u> 2 



Its finite temperature 



determinant is A(x m ,f3) = 2nhsin([3hui m )/Muj m . This 
goes through zero at [3 = ir/huj m and becomes negative 
for (3 > ir/huj m , thus signaling that x(t) = x m becomes 
unstable. At the same time, two new classical paths ap- 
pear, symmetric with respect to x m , as depicted in Fig. [l|. 
Therefore, at xq — x m , we go from a single-path regime 
to a triple-path regime as we cross (3 = Tr/fko m . The two 
new paths are degenerate minima, whereas the path that 
sits still at x m becomes a saddle-point of the action, with 
a single negative mode. 

As (3 grows beyond n/huj mi an analogous situation oc- 
curs for other values of xo inside the well. When the 
fluctuation determinant around the classical path for a 
given Xo vanishes, a new classical path appears. Its single 
turning point lies opposite, with respect to x mi to that of 
the formerly unique path. It may be interpreted as a tun- 
neling trajectory, since it traverses a classically forbidden 
region of V . If (3 is increased further, this tunneling tra- 
jectory splits into two, as illustrated in Fig. ||. One is 
a local minimum of the action, whereas the other is a 
saddle-point, with only one negative mode. Again, we 
have transitioned from a single to a triple-path regime. 
As j3 grows, the triple-path region spreads out around 



x m . The frontiers of that region are defined by the points 
x such that [<9t(x , y)/dy] y=x ^ = 0. 

The phenomenon just described is an example of catas- 
trophe. It takes place whenever the number of classical 
trajectories changes, with one or more of them becoming 
unstable, as we lower the temperature. Conversely, we 
may say that the phenomenon is characterized by the co- 
alescence of two or more of the classical trajectories, as 
we increase the temperature. This is akin to the ocur- 
rence of caustics in Optics || , where light rays play the 
role of classical trajectories and the action is replaced 
with the optical distance. 

If we denote by s\ the coordinate associated with the 
direction of instability of our example, the action can be 
viewed as S = S{s\, . . .\xq, f3). If we use the basis of 
eigenfunctions of the fluctuation operator around classi- 
cal paths, si will correspond to the direction along the 
one eigenfunction whose eigenvalue goes through zero, 
with the dots referring to all others. Catastrophe The- 
ory allows us to write down the "normal form" , Sn, 
of this generating function in the three-dimensional sub- 
space made up by the unstable direction of state space 
(the set of paths) and the two control variables (related 
to [3 and xq); it is this subspace which is relevant for the 
study of the onset of instability. As catastrophes are clas- 
sified by their codimension in control space, which here is 
two-dimensional, only those of codimension 1 (the fold) 
or 2 (the cusp) can occur. The pattern of extrema then 
leads to the cusp, whose generating function is 



5jv(si;ui,«i) 



uisf 



ViSi 
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where u\ and v±, the control parameters, are related to 
(3 and Xq, respectively. This catastrophe is defined by 



dS 
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d 2 S 
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dsi ds\ 







(9) 



Eliminating s\ from these equations, we can draw the 
bifurcation set in control space, as well as the pattern of 
extrema of Sn (Fig. ||). 

We can also plot the classical action (i.e., the action 
for classical trajectories) as a function of Xo for differ- 
ent values of /?, by exploiting the relation between the 
cusp and the swallowtail catastrophes. The latter has a 
generating function whose normal form is 



W(si; a, b, c) = sf + asf + bsf + csi 



(10) 



a, 6, c are control variables. The extremum condition 
is, then, dW/dsi = 0. The identifications Sn = —c/5, 
i*i = 3a and v\ = 26, will make the additional condition 
that defines the swallowtail, d 2 W/ds\ = 0, coincide with 
the requirement OSn /ds\ = 0. In the usual jargon, the 
bifurcation set for the swallowtail coincides with the equi- 
librium hypersurface of the cusp. We can, then, draw the 
graphs for Sn versus v±, for various values of u±, shown 
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in Fig. 0(a,b). Note that v\ is related to (x Q — x m ) and 
U\ to (f3 — ir/huj m ); for small values of both, we expect 
the relations to be linear. 

New tunnelings, accompanied by new catastrophes, 
will occur as we keep increasing (3. From (3 = 2ir/huj m , 
we start having tunneling trajectories with i(0) = 0, i.e., 
with one full period. For these, the determinant of fluctu- 
ations vanishes because the first, not the second bracket 
of (||) goes trough zero. However, the fluctuation op- 
erator around such tunneling trajectories already has a 
negative eigenmode. This follows from the fact that the 
zero-mode, given by i c i(r), has a node. Thus, it can- 
not be the ground-state for the associated Schrodinger 
problem. This means that a new catastrophe takes place 
along a second direction in function space. 

This new catastrophe is also a cusp. At x m , as (3 be- 
comes larger than 2Tr/hu) m , the solution that sits still ac- 
quires a second negative eigenvalue — thus becoming un- 
stable along a new direction in function space — and two 
other solutions appear. They both have a negative eigen- 
value along the first direction of instability and a posi- 
tive eigenvalue along the new direction of instability. If 
we combine the information from both directions, we find 
that we go from two minima and a saddle-point with only 
one negative eigenvalue (hereafter, a one-saddle) to two 
minima, two one-saddles and one two-saddle. The two 
one-saddles are degenerate in action and represent time- 
reversed periodic paths. As (3 increases beyond 2tt /huj m , 
the same phenomenon takes place for xq around x m : the 
one-saddles that already existed for (3 < 2ir /hu m in the 
three-path region around x m become unstable along a 
second direction, and a five-path region grows around 
x m , with two minima, two (time-reversed, degenerate in 
action) one-saddles and a two-saddle. 

The new cusp can be cast into normal form using the 
second direction of instability: 

Sn(s2] u 2 , v 2 = 0) = 4 + u 2 s 2 2 ~ (s\ + y-) 2 ■ ( n ) 

The absence of the linear term in Eq. ( |TT[ ) [compare with 
Eq. (||)] reflects the degeneracy in action of the two one- 
saddles. However, there is another degeneracy: since V 
does not depend on if x c \(t) is a solution of the Euler- 
Lagrange equations, so is x c i(t + tq). If x c \(t) is periodic, 
Xc\{t-\-tq) describes the same path — and so has the same 
action — , but with another starting point. This can be 
represented in Eq. ( |TT|) by the choice u 2 oc [(xq — x m ) 2 + 
k(2n/huj m - (3)}, with k > 0. See Fig. |(c). 

As we approach (3 — 3n /hto rn , the situation becomes 
similar to the one near (3 = tt /hw m . The difference is that 
we now have to deal with a tunneling path with more 
than one and less than two periods. Near j3 = An/huj m , 
two-period tunnelings intervene, a situation similar to 
that at f3 = 2n/?uv rn . The pattern which develops is 
depicted in Fig. bl 

Despite the fact that we keep adding new extrema as 



we lower the temperature, many of which represent tun- 
neling solutions, only two of these are stable (i.e., min- 
ima). In a semiclassical approximation in euclidean time, 
these are the only extrema we have to sum over, mean- 
ing that N(xq, (3) never exceeds two. The local minimum 
is a tunneling solution, whereas the global one becomes 
the unique solution in the high-temperature regime. For 
(3 < n/fkUm, there will be regions with either one or two 
minima of the classical action. The transition from a 
single-minimum to a double-minimum region occurs at 
values of xo where the fluctuation determinant vanishes 
due to the appearance of a caustic, thus leading to a sin- 
gularity in the integrand of Eq. (|^) . This is not a disaster, 
however, as this singularity is integrable. (Such a singu- 
larity is an artifact of the semiclassical approximation. 
It disappears in a more refined approximation fo]fl|l0f , 
in which one includes cubic and quartic fluctuations in 
the direction of function space where the instability sets 
in.) Explicit calculations for a number of potentials, and 
the derivation of Eq. (^), will be presented in a future 
publication. 
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FIG. 1. (a) Stable paths at x m for (3 < Tv/hiOm (1) and 
(3 > ir/hiOm (2 and 3); (b) Sketch of how the extrema change 
along the unstable direction in function space. 



FIG. 2. (a) Stable paths at x m for /3 < f3 c (1), f3 = (3 C (1 
and 2&3) and (3 > (3 C (1,2 and 3), Tr/huj m < f3 c < 2Tv/huj m ; (b) 
Sketch of how the extrema change along the unstable direction 
in function space. 

FIG. 3. Bifurcation set for the cusp; pattern of extrema 
shown schematically. 

FIG. 4. Evolution of the classical action as f3 changes, 
(a) < (3 < ir/hio m ; (b) Tv/hio m < (3 < 2-K/TiLo m ; (c) 
2n/huJ m < f3 < 3iv/hLU m - 



FIG. 5. Partition of control space into p-solution regions 
(p=l, 3,5,7,...). 
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